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Abstract 

To solve the relativistic bound-state problem one needs to systematically and simultaneously 
decouple the high-energy from the low-energy modes and the many-body from the few-particle 
states using a consistent renormalization scheme. In a recent paper we have shown that one such 
approach can be a combination of the coupled cluster method as used in many-body theory and the 
Wilsonian exact renormalization group. Even though the method is intrinsically non-perturbative, 
one can easily implement a loop expansion within it. In this letter we provide further support for 
this aspect of our formalism by obtaining results for the two-loop renormalized <^> 4 theory. We show 
that the non-unitary representation inherent in our method leads to an economic computation and 
does not produce any non-hermiticity in the relevant terms. 
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The Tamm-Dancoff approximation 1] was developed in the 1950's to describe a rela- 
tivistic bound state in terms of a small number of particles. It was soon revealed that the 
Tamm-Dancoff truncation gives rise to a new class of non-perturbative divergences, since the 
truncation does not allow us to take into account all diagrams at a given order in perturba- 
tion theory. On the other hand, any naive renormalization violates Poincare symmetry and 
the cluster decomposition property. Two very different remedies for this issue are the use of 
ight-front dynamics (see, e.g. |2j) and the application of the coupled cluster method (CCM) 
3(. However, both methods are too complicated to attack the issues in a self-consistent way. 
More recently, Glazek and Wilson 4] and independently Wegner p] introduced an elegant 
approach for such problems, the so-called similarity renormalization group (SRG). There 
are several problems with this approach: The Hamiltonian can not be systematically diag- 
onalized in particle number space, the computation is extremely complex and there is no 
efficient non-perturbative scheme. 
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Our method resembles the SRG approach, since we employ a similarity transformation 
to decouple the high-energy from the low-energy modes, leading to partial diagonalization 
of the Hamiltonian. Since we apply a double similarity transformation using linked-cluster 
amplitudes, we produce a biorthogonal representation of the problem, which is not neces- 
sarily unitary. There is a long tradition of such approaches in nuclear many-body theory 
[lfll | . One of our goals in this letter is to show that the non-hermit icity is not a problem, and 
leads to an economic computation. The method is intrinsically non-perturbative (it can be 
conceived as a topological expansion in the number of correlated excitations), but one can 
implement a loop expansion within it, and we shall concentrate on that aspect in this letter. 
By construction, we design a transformation which does not produce any fake divergences 
due to "small energy denominators" which plague old-fashioned perturbation theory in the 
Hamiltonian approach 11]. The Poincare invariance and the cluster decomposition prop- 
erty can be maintained in principle at every level of truncation regardless of regularization 
scheme, by requiring a set of decoupling conditions 

In this letter we supply further evidence for the efficiency of our method by computing 
the 2-loop renormalization of 4 theory. We only include a succinct discussion of the key 
elements of the mathematical framework; details can be found in Ref. |6j. Notice that our 
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formulation does not depend on the form of dynamics, i.e. on the choice of quantization 
(hyper-)plane. 

Consider a system described by a Hamiltonian H(A) which has, at the outset, a large 
cut-off A. We assume that generally the renormalized Hamiltonian H eG (X) up to scale A can 
be expressed as 

H cS (\) = H(\) + H c (\), (1) 

where Hc(X) is a "counterterm" . Our aim is to construct the renormalized Hamilto- 
nian by obtaining this counterterm. We define two subspaces, the model-space P : 
{\L) (g) \0,b) h ,L < /j,} and the complement-space Q : {\L) (g) (\H) - \0,b) h ) ,fi < H < A}. 
The ket |0, b)h is the bare high energy vacuum (the ground state of the free high-momentum 
Hamiltonian). The P-space contains interacting low-energy states and Q-space contains the 
orthogonal complement (the symbols \L) and \H) denote generic low- and high-energy states 
respectively). Our renormalization approach is based on decoupling of the complement space 
Q from the model space P by using a non-unitary transformation. The transformation of 
H(A) is defined by 

H = e s>,A) e -s(M)#( A ) e s(M,A) e -s'(M) = H{n) + 5H(ji, A), (2) 

where the operator S(S') is a functional of certain mapping operators between P- and 
Q-space |6[. By means of this transformation one may identify the effective interaction 
5H(fj,, A) containing the physics between the scale A and \i. According to our prescription 
we expand S(S') in terms of independent coupled cluster excitation /, 

S = S m ^— j , S m = -hCj, 

m=0 / 

^' = E^(x) m ' s' m =Y^ g i c i- ( 3 ) 

m=0 1 

Here the primed sum means that 1^0, and momentum conservation is included in si and 
s'j. The d and C\ are annihilation and creation operators in the high-energy Fock space 
for a given quantization scheme (e.g., equal time or light-cone). The indices / define a 
subsystem, or cluster within the full system of a given configuration. Notice, that the choice 
of the operators S(S') is not generally unique, due to Haag's theorem This ambiguity 
corresponds to the possibility of the choice of a different but equivalent representation of 
the canonical variables. However the physical quantities remain invariant under change of 
operators S(S'). 



It is well-known in many-body applications that the parametrization Eq. (jSJ) has the 
following properties: 1) it satisfies proper size-extensivity and conforms with the Goldstone 
linked-cluster theorem at any level of approximation 2) it is compatible with Hellmann- 
Feynman theorem, 3) the phase space {sj, s'j} for a given m is a symplectic different iable 
manifold. These features are related to each other and one can not give up one without losing 
the others as well Q]. The price we pay for these desirable features is that the representation 
is not longer hermitian. 

One may now impose the decoupling conditions, leading to diagonalization of the trans- 
formed Hamiltonian matrix, 

QHP = -f (0\C ie - § He § \0} = 0, (4) 
PHQ = -> (0\e § 'e- § He § e- § 'C\\0) = 0, (5) 

where I 4 0. One can show that the right hand side of Eq. (14151) are derivable from the 

n 

dynamics of the quantum system [6j. Thus the effective low-energy Hamiltonian is 

H cS = PHP = h (b, 0\e S '^ A) e- § ^ A) H(A)e S ^' A) e- § '^ A) \0, b) h . (6) 

One may prove that the effective low-energy operators Eq. (JBJ) supplemented by decoupling 
conditions Eq. (|4I5|1 indeed have the same low-energy eigenvalues as the original Hamiltonian 
The states in the full Hilbert space are constructed by adding clusters of high-energy 
correlation to states in the P-space, or equivalently by integrating out the high-energy modes 
from the Hamiltonian. It is immediately clear that states in interacting Hilbert space are 
normalized, due to the nature of a similarity transformation. 

The decoupling property makes the P sector of the truncated Fock space independent 
of the rest. This means that the contribution of the excluded sector is taken into account 
by imposing the decoupling conditions. One can show that the energy- dependent Bloch- 
Feshbach formalism [12| is thereby made free of the small-energy denominators which spoil 
perturbation theory 6(. One can then determine the counterterm by requiring coupling 
coherence ^| , namely that the transformed Hamiltonian Eq. (JHJ) has the generic form given 
in Eq. (Q), with A replaced by //. This requirement must be satisfied on an infinitely long 



1 It is of interest that according to our approach, various effective low-energy Hamiltonians can be con- 
structed without invoking perturbation theory or hermiticity. 
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renormalization group trajectory and thus produces a renormalized Hamiltonian. The indi- 
vidual amplitudes for a given m, {s™, s'f 1 } = {sj, s'j} m , have to be fixed by the dynamics of a 
quantum system incorporated the decoupling conditions Eq. (|4I5|) . Equations (|4I5|) provide 
two sets of exact, microscopic, operatorial coupled non-linear equations for the ket and bra 
states which describes the flow of the coefficients in Eq. (J3J). One can solve the coupled equa- 
tions in Eq. (0J) to obtain {s/} m and then use them as an input in (J5J). The conditions given 
in Eq. (@J) and (0) imply that all interactions of high-momentum particles should be removed 
from the transformed Hamiltonian H eS (fi) in Eq. ©. These are sufficient requirements to 
ensure partial diagonalization of the Hamiltonian in particle and momentum space. 

In practice one needs to truncate both sets of coefficients {§i, s'j} m at a given order of m 
in the ratio of cutoffs. A consistent truncation scheme is the so-called SUB(n, m) scheme, 
where we include up to n-body operators {S, S'} and truncate the expansion in (yu/A) at 
order m. The choice of n depends on the Hamiltonian interaction and needs to be fixed 
from the outset. 

We now apply this formalism to the computation of the effective Hamiltonian for <f) A 
iheory up to two- loop order in equal-time quantization. The bare 4 theory Hamiltonian is 

H = j (fx ()f 2 (x) + ( - V 2 + m 2 )0(x) + g<p\x^j . (7) 

According to our logic the ultraviolet-finite Hamiltonian is obtained by introducing coun- 
terterms, which depend on the UV cutoff A and some arbitrary renormalization scale. This 
redefines the parameters of the theory and defines the effective low-energy Hamiltonian. The 
renormalized Hamiltonian has the form 



H = j d z x (yvr 2 (x) + ~y/Zi</>(x)( - V 2 + Z m m 2 ) Jzl<\>{x) + Z^g^ix) + ... 



(8) 



Even though t 
is well known 



re newly generated interactions are sensitive to the regularization scheme (as 



16^ . a sharp cutoff may lead to new non-local interaction terms), nevertheless 
one can ignore these if they are irrelevant in the renormalization group sense. We now 
split field operators into high- and low-momentum modes; <ft(x) = <Pl{x) + 4>h{x), where 
denotes modes of low-frequency with momentum k < /i and 4>h(x) denotes modes 
of high-frequency with momentum constrained to a shell fi < k < A. The field 0_ff(x) can 
be conceived as a background to which the 0x,(x)-modes are coupled. Therefore, in the 
standard diagrammatic language, integrating out the high-frequency modes (Ph(%) implies 



that only high-frequency modes appear in internal lines. The field (Ph{x) is represented in 
Fock space as 

M*)= E -7^Ke^ + 4e-^), (9) 

where uj k = Vk 2 + m 2 and the operators a fc and a\ satisfy the standard boson commutation 
rules. From now on all summations are implicitly over the high-frequency modes // < k < A. 
The Hamiltonian in terms of high- and low-frequency modes can be written as, after normal 
ordering with respect to high-frequency modes, 

H = H 1 + H 2 + V B + V C + V A , (10) 

where we define, 

Hi = J Q7rI( a ;) + i0 L ( a ;)(-V 2 + m 2 )0 i (x)+^ 4 i (^, 
H 2 = y^fcOfeafc, 

/ e i(p+q+r-k)x Q e i(p+q-r—h)x 
—==a\a p a q a r + ——==a\a\a q a r 

e i(p+q-k)x J e i(k-p)x 

2a; r 

4444 + ^ 444 + ^ 44 + v c 4> 



V^B = 4 



V c = 


^Z 1 


V A = 




V c = 


/6(j> L {x) 
\ u p 


V c = 


2</> L (x)- 



ZU3 r / y/LUhLUp 



—ikx , i , „—i(k+p)x 



= ^ ; — . (11) 



The high-energy configurations in the Fock space are specified by {C/ — > rL=i a fc;} and 
{C] — > rii=i °L)}- Up to two-loop expansion, our renormalization scheme requires to keep 
S(S') at least to order n = 4, which allows us to eliminate the pure terms Vc and Va at 
a lower level of expansion. The S(S) operators consistent with a SUB(A,m) truncation 
scheme are, 

s m = Jj2 ( §1 4 + si 44 + s ™ 444 + ^ 4444) , 

S' m = J ^2 flfc + ^™ ° fea P + ^™ a fe«p4 + a k a p a q a r ^ . (12) 



One can expand Eqs. ()4|5|) in terms of fi/A, leading to the introduction of a consistent 
hierarchy of equations in powers of m p, which can be solved for the coefficients S m (S' m ) 
in Eq. ([12)1 . We split the diagonalization of the Hamiltonian matrix in an upper and lower 
triangle part, using the double similarity transformation. One may notice that the "most 
non-diagonal" terms in the Hamiltonian are Vc and Va (in the light-front Hamiltonian 
such terms do not exist because modes with longitudinal momentum identically zero are 
not allowed). The potential Vb is already partially diagonalized and does not change the 
vacuum of the high-energy states. Therefore, here we employ a minimal scheme, aiming at 
removal of Va and Vc only. 

We restrict ourselves to the elimination of the high-energy degrees of freedom up to the 
first order in the coupling constant g and second order in the ratio of cutoffs fi/A. Therefore, 
our truncation scheme is called SUB (4, 2). For m = one finds, 

V 1 V 2 
9 1 - 9 2 - -a- c 



UJk LUj- + UJ p ' 

V 3 V 4 

S 3 = -9 ^ , S 4 = -g ^ , (13) 

UJ k + UJ p + UJ q ' L0 k + UJ p + L0 q + UJ r 



where the are defied in Eq. Here, one has * = S< Q. At this s ta g e the results 

for the one-loop renormalization can be computed. We evaluate the effective Hamiltonian by 
substituting S(S') from Eqs. (fT2"j) and (|T3*j) into Eq. (JHJ). In order to achieve renormalization, 
one should identify the potentially divergent terms ( when A — > oo) in the expansion of 
H cS (fi). Such a process generally can be done by inventing a power-counting rule, using the 
property S n ~ ^S^-i- Here we take uo k ~ \k\ for /x ^> m and replace by J T^p-- The 
standard tadpole one-loop mass renormalization arises from Vb due to normal-ordering. We 
add this divergent term to Hi and renormalize the bare mass 

^l-loop = (0|Vs|0> = 6g J2j = ^( A 2 _ ^ j d S X(j) 2 (x)j 

1-^(A 2 -^ 2 )- (14) 



J rn 



In this order the contribution of the terms [Vc, 5], [Va, S'} and [Hi,S(S')] are zero, after 
projection on to the high-energy vacuum. The only divergent contributions come from 
[V^ 3 ' ) , S'o^ 3 ''] due to a double and third contraction of the high-frequency fields respectively. 
There are two other divergent terms, ([V^ , Sq^], however they are harmless and are 
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cancelled out by the divergence of [[H 2 , So], S' 2 ^]. One thus obtains, 

sh = 1892 f fMfM c i(k+p)(x-v) 

(2vr) 6 J u k u p (u k + u p ) 

_ f H^(v) ct{k+P+q)(x - y) _ (15) 

(27r) 9 J U k U p U q {u k + U p + U q ) 

In general evaluation of integrals like Eq. f)15|) may produce non-localities. This is due to the 
fact that the total momentum in integrands of Eq. (|15|h namely r\ = p + q and r 2 = k + p + q 
are in the low-momentum space. To evaluate such integrations, one can firstly reduce the 
potential divergent integrals by a change of variable, for example for the first integrand 
we use p,q —> p,r\, and then expand the integrand in r\jp. Therefore, after expansion 
and evaluating the momentum integrals, one may be faced with non-analytic terms in the 
low-momentum space. However here these are irrelevant and will thus be ignored. We find 

+ Y6^I ln (^) / rf3 a:(V0(x)) 2 + finite terms. (16) 
One can immediately deduce the renormalization factors Z g and from above expression 

The unknown coefficients in expression Si is computed by making use of Eq. (|13[) and solving 
coupled equations (HJ), therefore one may yield, 



Si = %= (20 L (x) - 2,vr L (x)0i(x) - W 1 , 
u^2u k \ Up ) u k ^ v\ 

S\ = , ^f^ P)X (1 _ i 2n L (x)4>M) - fa, Si] + ^VXS? 

(Uk + U n )\ U k U n \ J Uk+U n \ ' v\ 



sl = -7 9 rz , =7r L (x) - - 9 - , (v\st + V[V£, st~ l ] 

[U k + U p + U q Y^2u k UpU q [U k + Up + U q )\ ^ 



St = —, 9 - K Y\Vt v ,Sl]. (19) 



(U k + Up + U q + U 



f I 1 



In the above expression summation over dummy momentum indices is assumed. One can 
find S[ in the same manner by exploiting Eq. (JSJ and using Eq. (|T9*|) as an input, which 
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leads to 



with the notations, 



v= 1,...,4, 



(20) 



/la ^ 9_ 



/2a 



Si" 
S'i 

Cjf/3a 
g/4a 



l/=l ' u=l 

- ( E ^r 2)a ^( g ) - £ Iyr a 5? + [vi, sj- 

p l/=l ' I/=l 



jp t ujq 
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r Yin, s[ {4 ^ )a } 

{u k + u p + u q + u r ) ^ 



(21) 



The only divergent contribution up to order g 2 arises from, 



6H=-(0\[H 1 ,S 1 ],S' 



After the evaluation of the leading divergent part, we find that 



(22) 



5H 



392 i,/ A 



d 3 x7i 2 (x), 



(23) 



16vr 4 

which contributes to the two-loop wave-function renormalization Z n . By comparing Eqs. ()18|) 
and (|23|h one may conclude that Z n = ZT , as it should be. To finish the renormalization 
up to two-loop order, one should also take into account the contribution at order g 3 . The 
divergent terms at this level originate from 



5H = -(0\ [{V A + l/2Vc + V B ),S ],S' |0>. 



(24) 



After a straightforward but lengthy computation one can obtain the leading divergent parts, 



6H 



Til 
2vr 4 L 



d 3 x<j) 4 (x), 



(25) 



this term should be added to Eq. (JT3J) , therefore one can immediately deduce the correct 
total renormalization factor Z g up to two-loop order, 



1 + 



2tt 2 ~ \{i 



],!(-) 



,'A 
51 In ( - 



(26) 
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One can now immediately obtain the well-known [15| two-loop /3-function and anomalous 
dimension by making use of Eqs. (JT81 Efijl . 

It is important to point out that the diagonalization at first order in the coupling constant 
defines a correct low-energy effective Hamiltonian which is valid up to order g 3 . Having said 
that, from Eq. (J2U)) one can observe that the non-hermiticity of the S operator appears at 
order g 2 and in a lower order of fi/A. As we have shown, non-hermiticity is negligible up to 
two-loop order (asymmetric terms appear in irrelevant contributions). We conjecture that, 
for the present model, non-hermitian terms only appear in irrelevant contributions, whatever 
the order of truncation. 

One should note that Eqs. (JH |HJ) are not fully consistent with the Hellmann-Feynman the- 
orem |6j], although this choice considerably simplifies the diagonalization of the Hamiltonian 
operator. This potentially leads to a different truncation scheme for the renormalization of 
other operators. The decoupling requirements in the exact form Eq. ()4I5|) preserve Poincare 
symmetry and the cluster decomposition property regardless of regularization used jf|. Since 
this is a continuous symmetry, its preservation leads to an infinite set of constraints on the 
phase space, which has been coded in the decoupling equations. Hence, it is of interest to 
consider the sensitivity of the Poincare symmetry with respect to a given truncation scheme 
in the light-front dynamics. 

In this letter we have employed a sharp cutoff, however this idealization should be removed 
since generally it may lead to pathologies in renormalization, since it induces non-locality and 
moreover potentially violates the gauge symmetry. We have only investigated a perturbative 
approach to the problem. Our method is in principle non-perturbative; this aspect remains 
to be exploited. 
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and UMIST grant. The work of NRW is supported by the UK engineering and physical 
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